Two-particle correlations on transverse momentum and momentum dissipation in 

Au-Au collisions at ^snn = 130 GeV 
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Measurements of two-particle correlations on transverse momentum pt for Au-Au collisions at 



^SNN ~ 130 GeV are presented. Significant large-momentum-seale correlations are observed for 
charged primary hadrons with 0.15 < pt < 2 GeV/c and pseudorapidity \ri\ < 1.3. Such correlations 
were not observed in a similar study at lower energy and are not predicted by theoretical collision 
models. Their direct relation to meam-pt fluctuations measured in the same angular acceptance is 
demonstrated. Positive correlations are observed for pairs of particles which have large pt values 
while negative correlations occur for pairs in which one particle has large pt and the other has much 
lower Pt- The correlation amplitudes per final state particle increase with collision centrality. The 
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observed correlations are consistent with a scenario in which the transverse momentum of hadrons 
associated with initial-stage semi-hard parton scattering is dissipated by the medium to lower pt. 

PACS numbers: 25.75.-q, 25.75.Gz 



I. INTRODUCTION 

Studying two-particle correlations and event-wise fluc- 
tuations can provide essential information about the 
medium produced in ultrarelativistic heavy ion colli- 
sions [U, [3, 0] ■ At the collision energies available at the 
Relativistic Heavy Ion Collider (RHIC) energetic parton 
scattering occurs at sufiicient rate to enable quantita- 
tive studies of in-medium modification of parton scatter- 
ing and the distribution of correlated charged hadrons 
associated with those energetic partons. Modification 
of those correlation structures is expected as the bulk 
medium produced in ultrarelativistic heavy ion collisions 
increases in spatial extent and energy density with in- 
creasing collision centrality. Analyses of the central- 
ity dependence in Au-Au collisions of high-pt back-to- 
back jet angular correlations based on a leading-particle 
technique {e.g., leading-particle pt > 4 GeV/c, associ- 
ated particle pt < 4 GeV/c) reveal strong suppression 
for central collisions [3, suggesting the development 
of a medium which dramatically dissipates momentum. 
Complementary studies of the lower-pt bulk medium, its 
correlation structure on transverse momentum, and how 
those correlations evolve with collision centrality pro- 
vide a measure of the momentum transport from the few 
GeV/c range to lower pt of order a few tenths of a GeV/c 
where the bulk hadronic production occurs. Such studies 
are an essential part of understanding the nature of the 
medium produced in heavy ion collisions at RHIC. 

In addition to jet angular correlations at high-pt sub- 
stantial nonstatistical fluctuations in event-wise mean 
transverse momentum (pt) of charged particles from 
Au-Au collisions were reported by the STAR 6] and 
PHENIX [Tj] experiments at RHIC. {pt) fluctuations at 
RHIC are much larger than those reported at the CERN 
Super Proton Synchrotron (SPS) with one-tenth the CM 
energy J8| , and were not predicted by theoretical mod- 
els [y, S [13, EH- {Ptj fiuctuations could result from 
several sources including collective fiow {e.g., elliptic 
flow [1^ when azimuthal acceptance is incomplete), lo- 
cal temperature fluctuations, quantum interference [l3|, 
final-state interactions, resonance decays, longitudinal 
fragmentation |14| . and initial-state multiple scatter- 



ing [15| including hard parton scattering ^9] with sub- 
sequent in-medium dissipation [l^. {pt) fluctuations can 
be directly related to integrals of two-particle correlations 
over the pt acceptance. Correlations on pt, by provid- 
ing differential information, better reveal the underlying 
dynamics for the observed nonstatistical fluctuations in 

In this paper we report the first measurements at RHIC 
of two-particle correlations (based on number of pairs) 
on two-dimensional (2D) transverse momentum space 



{pti,Pt2) for all charged particles with 0.15 < pt < "2 
GeV/c and \ri\ < 1.3 (pseudorapidity) using the ^s^n = 
130 GeV Au-Au colfisions observed with the STAR detec- 
tor [13 . This analysis is intended to reveal the response 
of the bulk medium to strong momentum dissipation and 
probe the dynamical origins of (pt) fiuctuations. The 
data used in this analysis are described in Sec. [TT] and 
the analysis method, corrections and errors are discussed 
in Sec. IIIII Models and fits to the data are presented in 
Sees. llVlandlVl respectively. A discussion and summary 
are presented in the last two sections I VI I and IVIII 



II. DATA 

Data for this analysis were obtained with the STAR de- 
tector |17| using a 0.25 T uniform magnetic field parallel 
to the beam axis. A minimum-bias event sample (123k 
triggered events) required coincidence of two Zero-Degree 
Calorimeters (ZDC); a 0-15% of total cross section event 
sample (217k triggered events) was defined by a thresh- 
old on the Central Trigger Barrel (CTB) scintillators, 
with ZDC coincidence. Event triggering and charged- 
particle measurements with the Time Projection Cham- 
ber (TPC) are described in [l3|. Approximately 300k 
events were selected for use in this analysis. A primary 
event vertex within 75 cm of the axial center of the TPC 
was required. Valid TPC tracks fell within the detector 
acceptance used here, defined by 0.15 < pt < 2.0 GeV/c, 
\r]\ < 1.3 and 27r in azimuth. Primary tracks were defined 
as having a distance of closest approach less than 3 cm 
from the reconstructed primary vertex which included a 
large fraction of true primary hadrons plus approximately 
7% background contamination from weak decays and 
interactions with the detector material. In addition ac- 
cepted particle tracks were required to include a mini- 
mum of 10 fitted points (the TPC contains 45 pad rows 
in each sector) and, to eliminate split tracks {i.e., one 
particle trajectory reconstructed as two or more tracks), 
the fraction of space points used in a track fit relative 
to the maximum number expected was required to be 
> 52%. Particle identification was not implemented but 
charge sign was determined. Further details associated 
with track definitions, efficiencies and quality cuts are 
described in 18, 19,]. 



III. DATA ANALYSIS 
A. Analysis Method 

Our eventual goal is to determine the complete struc- 
ture of the six-dimensional two-particle correlation for 



4 



all hadron pair charge combinations. Toward this goal 
the two-particle momentum space was projected onto 2D 
subspace {pti,Pt2) by integrating the pseudorapidity and 
azimuth coordinates (771, ?72, 0i, ^2) over the detector ac- 
ceptance for this analysis, \ri\ < 1.3 and full 2tt azimuth. 
Projection onto 2D subspace {pti,pt2) is achieved by fill- 
ing 2D binned histograms of the number of pairs of parti- 
cles for all values of 77, (j) within the acceptance. Comple- 
mentary correlation structures on relative pseudorapidity 
and azimuth coordinates with integration over transverse 
momentum acceptance are reported in [20I [2l| . 

The quantities obtained here are ratios of normalized 
histograms of sibling pairs (particles from the same event) 
to mixed-event pairs (each particle of the pair is from a 
different, but similar event) in an arbitrary 2D bin with 
indices a, h representing the values of pti and pt2 (see 
discussion below). The normalized pair- number ratio fab 
introduced in 22] is here defined by 

fah = f^ ab, sib / flab, mix ^ (1) 

where nab,sib = nab,sib/ Y,ab^ab,sib (sum over aU 2D 

bins), flab^raix — '^ab,rnix / ^^^b^*^^^^'^^ ^ and ?T'afe,sib and 

nab,mix are the inclusive number of sibling and mixed- 
event pairs, respectively, in 2D bin a, b. Histograms and 
ratios fab were constructed for each charge-sign combi- 
nation: (+,+), (—,—), (+,— ) and (—,-!-). Ratio fab is 
approximately 1, while difference {fab — 1) measures cor- 
relation amplitudes and is the quantity reported here. 

The exponential decrease in particle yield with in- 
creasing pt degrades the statistical accuracy of fab at 
larger transverse momentum, thus obscuring the statis- 
tically significant correlation structures there. In or- 
der to achieve approximately uniform statistical accuracy 
across the full pt domain considered here, nonuniform 
bin sizes on pt were used. This was done by noting that 
the charged hadron pt distribution, dN/ptdpt, for Au-Au 
collisions at ^/snn — 130 GeV is approximately expo- 
nential for 0.15 < Pt < 2 GeV/c flS'l and by dividing 
the running integral of that exponential distribution into 
equal bin sizes. This procedure provides a convenient 
mapping from pt to function X(pt) = 1 — exp{— (rn^ — 
mo)/0.4 GeV} where < X < 1, mt — \/pl+rrif^, and 
mo (here assumed to be the pion mass m^) is a map- 
ping parameter from coordinate pt to X [23]. Equal 
bin sizes in X therefore have approximately the same 
number of sibling pairs. For this analysis 25 equal 
width bins on X from X{pt = 0.15 GeV/c) = 0.15 to 
X(j,t = 2.0 GeV/c) = 0.99 were used 

Normalized pair-number ratios were formed from sub- 
sets of events with similar centrality (multiplicities dif- 
fer by < 100, except < 50 for the most-central event 
class) and primary-vertex location (within 7.5 cm along 
the beam axis) and combined as weighted (by sibling 
pair number) averages within each centrality class. The 
normalized pair-number ratios for each charge-sign were 

combined to form like-sign (LS: H — h, ) and unlike- 

sign (US: H — , — V) quantities. The final correlations re- 
ported here were averaged over all four charge-sign quan- 



tities, resulting in the correlation structures common to 
all charge-sign combinations. Hence we refer to these fi- 
nal results as charge-independent (CI = LS + US) corre- 
lations even though they are constructed from quantities 
which depend on the charge signs of the hadron pairs. 
The correlation measure reported here is therefore the 
CI combination for f\X{pi\),X{pt'2S\ — 1- 

Deviations of event- wise (pt) fluctuations from a 
central-limit-theorem reference [25l [26j are measured by 
scale-dependent [i.e., 77, (/) bin sizes) variance difference 
Ao-p^.„ introduced in [0], where it was evaluated at the 
STAR (?7, 0) detector acceptance scale. Act^^.^ can be ex- 
pressed as a weighted integral on (j)ti,Pt2) of pair-density 
difference psib — Pmix, where two-particle densities psib 
and Pmix are approximated by the event-averaged num- 
ber of sibling and mixed-event pairs per 2D bin, re- 
spectively. Both densities are normalized to the event- 
averaged total number of pairs. A(Tp^.„ can be rewritten 
exactly as a discrete sum over pt products |25|] [first line 
in Eq. ^ below], and the summations approximated in 
turn by the weighted integral of the pair density differ- 
ence [second line in Eq. ^] according to, 

11^ 

^'^If.n = ^-HY^ (ptnPtje-Pt) (2) 

j=l i^i' = l 

~ J J dptldpt2PtlPt2 {Psib - Pmix) 

= plN{r{ptuPt2)-l), 

where weighted average {r{pti,Pt2) ~ 1) is defined in the 
last line with weight ptiPt2PmixiPti,Pt2) and the integral 
of Pmix, J J dptidpt2PtiPt2Pmix is Pt [13]. In Eq._([2|) 
Nj is the event- wise number of accepted particles, N is 
the mean of Nj in the centrality bin, e is the number of 
events, j the event index, pt is the mean of the ensemble- 
average Pt distribution (all accepted particles from all 
events in a centrality bin), and i,i' are particle indices. 
Eq. ^ relates nonstatistical {pt) fiuctuations at the ac- 
ceptance scale to the weighted integral of psib — Pmix, the 
latter difference being related to the two-particle num- 
ber correlation density. In the present analysis we mea- 
sure normalized pair-ratio distributions f[X{pti), X{pt2)\ 
exhibiting two-particle number correlations on pt which 
correspond to excess {pt) fluctuations. 

B. Corrections and centrality 

Corrections were applied to ratio f for two-particle 
reconstruction inefficiencies due to overlapping space 
points in the TPC (two trajectories merged into one 
reconstructed track) and intersecting trajectories which 
cross paths within the TPC and are reconstructed as 
more than two tracks. These corrections were imple- 
mented using two-track proximity cuts (28j at various ra- 
dial positions in the TPC in both the longitudinal (drift) 
and transverse directions (approximately along the pad 



5 



rows). The track pair cuts were applied to both psib 
and Pmix as in HBT analyses [isj . Small- momentum- 
scale correlation structures due to quantum interfer- 
ence, Coulomb and strong final-state interactions [l^ 
were suppressed by eliminating sibling and mixed-event 
track pairs (~3% of total pairs) with I771 — 772] < 0.3, 
- 02| < 7r/6 (azimuth), \pti - Pt2\ < 0.15 GeV/c, 
if Pt < 0.8 GeV/c for either particle. The small- 
momentum-scale correlation (SSC) structures are most 
prominent in the lower-p^ domain of the 2D {pti,Pt2) 
space along the pti = Pt2 diagonal and were shown to 
be similar in amplitude and location to simulations (29j 
which account for quantum interference correlations and 
Coulomb final-state interaction effects usin g pa ir weights 
determined by HBT analyses for these data [13] . The pre- 
ceding cuts were optimized flO^ to eliminate SSC struc- 
ture without affecting the large-momentum-scale corre- 
lation (LSC) structure which is of primary interest here. 
The track-pair cuts generally have small effects on the 
LSC; uncertainties which result from application of these 
cuts are discussed in Sec. IIII CI and are negligible com- 
pared to the large momentum scale structures studied 
here. 

Four centrality classes labeled (a) - (d) for central to 
peripheral were defined by cuts on TPC track multiplicity 
A*' within the acceptance by (d) 0.03 < N/Nq < 0.21, (c) 
0.21 < iV/7Vo < 0.56, (b) 0.56 < N/Nq < 0.79 and (a) 
N/Nq > 0.79, corresponding respectively to approximate 
fraction of total cross section ranges 40-70%, 17-40%, 5- 
17% and 0-5%. Nq is the end-point ^3d\ of the minimum- 
bias multiplicity distribution. 

The centrality dependence of quantity f — 1 is shown in 
Fig. [T] as perspective views for the four centrality classes 
used here. This correlation measure represents the num- 
ber of correlated particle pairs per final-state pair in each 
2D bin, and therefore contains a dilution factor 1/N rel- 
ative to the LSC measure presented in [2l|, N{f — 1) 
whose amplitudes are of order one. The structures in 
Fig. [1] are therefore numerically a few permil for cen- 
tral Au-Au collisions but are highly significant statisti- 
cally as seen by comparing to the statistical errors. The 
dominant features in Fig. [1] are 1) a large-momentum- 
scale correlation 'saddle' structure with positive curva- 
ture along the X{pt)j: = X{pti) + X{pt2) sum direction 
from [X{pti), X{pt2)] (0,0) to (1,1) and a correspond- 
ing negative curvature along the X(pt)A = X{pti) — 
X{pt2) difference direction from [X (pti) , X {pt2)] (0,1) to 
(1,0), and 2) a narrow peak structure at large X{pt) 
{pt > 0.6 GeV/c). With increasing centrality the neg- 
ative curvature of the LSC saddle shape along the differ- 
ence variable increases in magnitude, the positive curva- 
ture along the sum variable decreases, and the magnitude 
of the peak at large X{pt) also decreases. Without the 
SSC cuts a relatively small peaked structure with am- 
plitude of order 0.004 (peripheral) to 0.0005 (central) 
is present for X{pt) < 0.3 {pt < 0.25 GeV/c) which 
weakens in amplitude but visibly persists to X{pt) < 0.6 
{pt < 0.5 GeV/c). 




FIG. 1: Symmetrized pair-density net ratios 

f[X{pti),X(pt2)] — 1 for all nonidentified charged pri- 
mary particles for (a) most-central, (b) mid-central, (c) 
mid-peripheral, and (d) peripheral Au-Au collision events at 
^/sNN = 130 GeV/c. Note the scale change for panels (c) 
and (d) and auxiliary pt scale in units GeV/c in panel (a). 
SSC were removed using track pair cuts (see text). Errors 
are discussed in Sec. IIII CI 



An upper limit estimate for resonance contributions 
was obtained using Monte Carlo simulations [1^ assum- 
ing 70% of the primary charged particle production is 
from resonance decays. The correlations were simulated 
by populating the events with a sufficient number of p'^, w 
two-body decays to account for 70% of the observed mul- 
tiplicity. These two-body decay processes produced a 
small saddle-shape correlation with curvature opposite 
to the data and amplitude at the corners approximately 
0.0002 for the most-central data, increasing as l/N for 
the remaining centrality bins. The saddle-shape struc- 
tures in Fig. [1] cannot be explained with resonance de- 
cays. 

The same analysis applied to Pb-Pb collisions in 1.1 < 
Ucm < 2.6 at the CERN SPS did not reveal any sta- 
tistically significant CI correlations [3l| when SSC (see 
Sec. IIII Bp were removed with pair cuts. The analysis in 
[32i] of proton -I- proton and various nucleus -I- nucleus 
collision data from the CERN SPS for 1.1 < y^n < 2.6 
without those pair cuts revealed SSC peaks at low X{pt) 
along the X{pt)^ direction. 
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C. Error analysis 

Per-bin statistical errors for f — 1 in Fig. [1] range 
from ^6-9% of the maximum correlation amplitude for 
each centrality [typically 0.00015, 0.00011, 0.00035 and 
0.001 for centralities (a)-(d) respectively] and are ap- 
proximately uniform, by design, over the 2D domain on 
X{pt). Statistical errors for N{r - 1) (~0.1 - 0.15) are 
less dependent on centrality. 

Systematic errors were estimated as in [2l| and 
are dominated by the 7% non-primary background con- 
tamination whose correlation with primary parti- 
cles is unknown. The upper limit on the systematic 
error from this source was estimated by assuming the 
number of correlated pairs associated with background- 
primary pairs of particles could range from zero up to 
the amount which would occur among 7% of the pri- 
mary particles and the remaining primaries. This con- 
servative assumption produces an overall ±7% uncer- 
tainty relative to the correlation amplitudes in Fig. [1] 
throughout the domain for X(pti,2) > 0.4. This error in- 
creases to ±16% at lower X{pt) where the contamination 
fraction is larger and is about ±12% in the off-diagonal 
corners of the [X (pti) , X {pt2)] domain. Multiplicative 
factors for quantity f — 1 which correct for the non- 
primary background contamination range from 1.0, as- 
suming background-primary particle pairs are correlated 
and increase both hab,sib and fiab.mix by 2 x 7% = 14%, 
to 1.14 if background-primary particle pairs are uncorre- 
lated but the non-primary background contributes 14% 
to nab,mix- Multiplication of the data in Fig. [T]by aver- 
age factor 1.07 provides an estimate of the background 
corrected correlation amplitudes. 

Additional sources of systematic error were evalu- 
ated. Uncertainty in the two-track inefhcency correc- 
tions have modest effects along the X{pti) = X{pt2) di- 
agonal (< 2%) and are negligible elsewhere. Tracking 
anomalies caused when particle trajectories intersect the 
TPC high-voltage central membrane significantly affect 
the X{pti,2) < 0.2 domain corresponding to the single bin 
at lowest X{pt), and the diagonal bins for X{pti^2) < 0.4 
by 20%. Final multiplicative correction estimates (not 
applied in Fig. [1]) and total systematic errors for f — 1 
varied respectively from 1.07 and ±7% for X{pti.2) > 0.4 
up to 1.16 and ±16-20% for X{pti,2) < 0.4 and 1.12 and 
±12% in the off-diagonal corners [i.e., near (0,1) and 
(1,0)]. 

Other potential sources of systematic error were stud- 
ied and determined to have negligible effects including 
primary vertex position uncertainty perpendicular to the 
beam direction, variation of tracking acceptance and effi- 
ciency with primary vertex location along the axis of the 
TPC, TPC drift speed and/or timing-offset fluctuation, 
sporadic outages of TPC read-out electronic components, 
angular resolution, multiplicity and primary vertex po- 
sition bin sizes used for producing mixed events, and 
charge sign dependence of the tracking efficiency. Con- 
version electron contamination is suppressed by the lower 



Pt acceptance cut and also by the pair cuts described in 
Sec. IIIIBI and also makes negligible contribution to the 
systematic error. 



IV. MODELING ONE- AND TWO-PARTICLE 
DISTRIBUTIONS ON pt 

Two features dominate the data in Fig. [TJ 1) a large- 
momentum scale saddle shape and 2) a peak at large 
X{pti) and X{pt2). In this section results from Monte 
Carlo collision models are analyzed in order to gain in- 
sight into the dynamical origin(s) of these two correlation 
structures in the data. Based on this study an analytical 
function is obtained which accurately describes the sad- 
dle shape and in Sec. |V] this function is used to fit the 
2D correlation data. 

The high energy nuclear collision model hijing 
which includes longitudinal color-string fragmentation 
and perturbative quantum chromodynamics (pQCD) 
based jet production and jet quenching, exhibits signifi- 
cant correlation structure at higher [X {pti) + X {pt2) > 
1.6] as shown in the left-hand panel of Fig. [5] for cen- 
tral Au-Au collisions. The predictions, which include jet 
production with jet quenching (default parameters) are 
qualitatively different from the data in Fig. [l] failing to 
produce any saddle-shape, but suggest the type of corre- 
lation structure produced by jets. The general structure 
of the HIJING predictions suggests that the peaks in the 
data at higher X{pt) are at least partly due to initial- 
state partonic scattering and fragmentation. Other the- 
oretical models which combine initial-state parton scat- 
tering, energy loss, dissipation, rescattering and recombi- 
nation 111, 3J] may eventually explain these correlation 
data, but relevant predictions are not available at this 
time. 
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FIG. 2: Symmetrized pair-density ratio r[X (pti) , X (prz)] — 
1 for unidentified charged particles and for central Au-Au 
collisions. Left panel: Default HiJiNG [9|] with jet quenching. 
Right panel: a Monte Carlo model [2^ which simulates event- 
wise global temperature fluctuations (see text). 



The saddle-shape correlation spans the entire momen- 
tum scale studied here, suggesting event-wise fluctuations 
of global event characteristics {e.g. temperature and/or 
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collective velocity of the bulk medium) as a possible 
source. If heavy ion collisions at RHIC thermalize then 
an ensemble of collision events would be characterized 
by a distribution of event-wise equilibrium temperatures 
reflecting event-to-event fluctuations in the initial condi- 
tions and time evolution of each colliding system. Based 
on this idea the transverse momentum correlations for an 
ensemble of such events can be predicted using a Monte 
Carlo model in which charged particle production is gen- 
erated by sampling the inclusive single-particle {pt,ri,4>) 
distribution obtained from the data. At mid-rapidity the 
inclusive distribution on pt for 0.15 < pt < 2 GeV/c 
is well approximated by exi>(—mt/T ) = exp(— /3mt) [18| 
where T is an effective temperature [3J] or inverse slope 
parameter and f3 — 1/T. Events were generated by 
sampling exp(— mt/T) where T fluctuates randomly from 
event-to-event according to a gaussian distribution about 
mean value Tq = l//3o; To was determined by the mea- 
sured Pt spectrum. 

The result of this Monte Carlo model for central Au- 
Au collisions at 130 GeV is shown in the right-hand panel 
of Fig. [2] where the mean and standard deviation (gaus- 
sian sigma) of the event-wise temperature distribution 
are Tq = 200 MeV and ctt/To = 1.5%. The predicted 
correlations are not sensitive to Tq but the overall cor- 
relation amplitude is directly sensitive to ut/Tq which 
was adjusted to approximate the overall amplitude of 
the data in Fig. [Ija). The global temperature fluctu- 
ation model accurately describes the saddle-shape. An 
analytical function based on this approach is derived in 
the remainder of this section and is used in the following 
section (Sec. |V| to fit the data. 

We seek an analytical representation of the LSC 
saddle-shape structure of the data in Fig. [1] that is 
both mathematically compact and physically motivated 
in order to conveniently characterize the centrality de- 
pendence and to infer thermodynamic properties of the 
medium. The above Monte Carlo results indicate that 
a successful representation should involve an averaging 
over the inverse slope parameter. In general the inverse 
temperature (3 can vary from event-to-event as well as 
internally within each event, reflecting the possibility of 
relative "hot spots" and "cold spots" in the final-state 
particle distributions. The number, location in source 
coordinates {e.g., r\z — space-time rapidity [s^l and 
— azimuth), amplitude, and angular extent of these per- 
turbations in /3 may vary for each event. In addition, 
for realistic collision systems both thermal and collective 
motions are involved such that parameter /3 becomes an 
inverse effective temperature [34] where fluctuations in /3 
could result from fluctuations in the local temperature 
of the flowing medium, the collective flow velocity itself, 
or a combination of both effects [s^. Event- wise effec- 
tive temperature is therefore represented by distribution 
T{rjz,^p) on source coordinates rjz and Lp, and similarly 
for I3{f]z,f)- 

The momentum of a particle at riz,ip in the final stage 
of the collision system is obtained by sampling thermal 



distribution eyip[—mt/T{r]z,ip)] = exp[—P{rjz,(p)mt] as 
illustrated in the diagram in Fig. EKa). In general the 
histogram of sampled T(rjz,(p) or P{r]z,ip) values for all 
particles in all events in the event ensemble, gi(/3), could 
be like the generic peaked distribution in Fig. [Sfb) with 
mean Pq and standard deviation a p. The inclusive nit 
distribution is then obtained by convoluting thermal dis- 
tribution exp[— /3(mt — Too)] with (?i(/3) given by 



dN 
mtdmt 



= A 



-0(mt-ma) 



(3) 



where A is a normalization constant. The global tem- 
perature fluctuation model is recovered when T{'r]z, f) is 
independent of source coordinate but varies from event- 
to-event. 

In the Monte Carlo model event-wise T — 1/(3 was ob- 
tained by sampling a gaussian distribution. It is therefore 
reasonable to represent gi(/3) by a peaked distribution 
which is here assumed to be a gamma distribution [s^ 
in order to obtain an analytic solution of the integral in 
Eq.© given by 



dN 
mtdmt 



A/[l -f /3o(mt - mo)/n/i„ct]"^'"" , (4) 



a Levy distribution [37], where l/ufiuct — '^fs/Po is the 
relative variance of gi{(3). The finite width of gi{(3) pro- 
duces a net increase in the yield at higher mj as illus- 
trated in Fig. [3{c). We emphasize that any finite- width 
peaked function gi{(3) results in an mt distribution which 
decreases less rapidly with increasing mt than thermal 



spectrum e 



The assumption of a gamma distri- 



bution for 5i(/3) is therefore not essential but is used for 
mathematical convenience and is justified by the capabil- 
ity of the mt distribution in Eq. ^ to describe the inclu- 
sive data. We note however that deviations of the mea- 
sured mt distribution from a thermal spectrum, quanti- 
fied by exponent n in the power-law mt distribution , 
can result from transverse expansion ^34 ] in addition to 
local and event-to-event fluctuations in /3(?7z, ^p) assumed 
in deriving Eqs. ([3]) and Consequently, fitting the 
1 /mt dN/ dmt spectra to obtain the power-law exponent 
n cannot by itself determine the relative variance of the 
effective temperature distribution, l/ufiuct, which is re- 
lated to the degree of equilibration. 

Similarly the two-particle distribution on {mti,mt2) is 
obtained by convoluting the two-particle thermal distri- 
bution exp[— /3i(TOti — mo)] exp[— /32(to(2 — tuq)] with the 
2D distribution of pairs of inverse effective temperature 
parameters (/5i,/32), where particles 1 and 2 sample lo- 
cal thermal distributions determined by P{rizi,ipi) and 
P{'ilz2, 'P2), respectively [see Fig. [31[a)]. The distribution 
of {Pi, (32) for all pairs of particles used in all events in 
the ensemble defines a 2D histogram and 2D distribution, 
92{f3i, Pi), illustrated in Fig. [3l panels (d)-(f) for three 
hypothetical cases. If the event ensemble distribution on 
P has finite width {up > 0), but are point-to-point uncor- 
related within each event, then (72 (/^i, P2) is symmetric on 
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FIG. 3: Diagrams illustrating the temperature fluctuation 
model. Panel (a): Source coordinates with two final state 
particles sampling local inverse temperatures /3i = (3{rjzi,ipi) 
and P2 = /3(??z2, <(52). Panel (b): Distribution gi{P) of sam- 
pled /3 values for all particles in all events of a centrality 
bin with mean f3o and standard deviation a p. Panel (c): 
Thermal model inclusive charged particle yield dN/mtdmt 
at mid-rapidity versus mt — mo with no temperature fluc- 
tuations {afj — 0, solid line) and with temperature fluctu- 
ations ((j/3 > 0, dashed line). Panel (d): 2D distribution, 
52 , /32 ) , of sampled pairs I3{r)zi,tpi) and P{riz2,'fi2) when 
there are no point-to-point temperature correlations within 
each source but large temperature variations within each 
event (non-equilibrium sources). Panel (e): Same as (d) ex- 
cept for global temperature fluctuations where each event is 
equilibrated but the equilibrium temperature fluctuates from 
event-to-event. Panel (f): Same as (d) except point-to-point 
temperature correlations occur within each event as evidenced 
by the positive covariance of distribution g2{Pi,p2)- 



Pi vs P2 (zero covariance) as shown in Fig.[3]Jd). For un- 
correlated /3 fluctuations or for mixed-event pairs, (?2 fac- 
torizes as g2iPi, P2) = <?i(/3i)3i(/32), implying zero covari- 
ance. On the other hand, if every event is thermally equi- 
librated, then each pair of particles from a given event 
samples the same value of (3 where f3i = (32- For this case 



(global temperature fluctuation model) (72(/3i,/?2) limits 
to a diagonal line distribution illustrated in Fig.EJe) and 
given by 52(/3i,/32) oc g[{Pi)S{Pi - (32), where 6{Pi - (32) 
is a Dirac delta-function. In this case g2 has maximum 
covariance and represents the conventional picture of an 
ensemble of equilibrated events with event-wise fluctu- 
ations in the global temperature. In general g2i(3i,(32) 
may have an intermediate covariance as illustrated in 
Fig. ^i)- In this case if (?2(/3i,/?2) is expressed as a 
product of a gamma distribution on the sum direction, 
/3e = /3i + 132 multiplied by a gaussian on /3a = f3i — P2 
(for mathematical convenience) , then an analytic expres- 
sion for the two-particle distribution results, given by a 
2D Levy distribution 



Fsib oc 1 + — 

2ns / 



1 - 



2nA + PanitT. 



(5) 



on sum and difference variables mtT. = friti + '7it2 ~ 2too 
and TTifA = TT^ti — 'mt2- Inverse exponents l/ris and 
I/tt-a are the relative variances of g2{(3i,f32) along sum 
and difference variables /3s and /3a respectively, and 
A(l/n)tot = l/ns — l/riA is the relative covariance 
of (72 [38l |. measuring velocity/temperature correlations. 
For the examples in panels (d), (e) and (f) of Fig. [31 
1/ns = 1/?T-A, 1/ns > and l/riA = 0, and > 
1/riA > 0, respectively. Mixed-event pair distribution 
Fmix{pti,Pt2), a product of one-dimensional Levy dis- 
tributions [Eq. (U))], has the form of Eq. ^ but with 

= UA = rifiuct- 
Ratio 



(6) 



referred to as a 2D Levy saddle, predicts a saddle-shape 
when 52(/3i,/32) has non-zero covariance and is the ana- 
lytical quantity to be compared to data. It can be tested 
by comparison to the data in Fig. [1] via chi-square fits. 
We emphasize for this 2D case that any peaked function 
52(/3i,/32) with nonzero covariance results in a 2D saddle 
shape distribution for rmodei- The gamma distribution 
times gaussian 2D model for 32 was chosen for mathemat- 
ical convenience but it is reasonable given the form of the 
measured event-wise (pt) distribution. The variance of 
g2 along the difference direction /3a measures the average 
degree of equilibration of the events in the ensemble. Rel- 
ative variance differences A(l/n)s = (l/^s ~ ^/nfiuct) 
and A(l/n)A = (I/^a — l/nfiuct) measure the saddle 
curvatures of r„iodei (and hence the data) along sum and 
difference directions at the origin, and are the quantities 
best determined by these fits. Sensitivity to the magni- 
tudes of the relative variances 1 /ns and 1 /?ia is discussed 
in the next section. 



V. ANALYTICAL MODEL FITS 

Data in Fig. [T] (excluding peak region X{pt)j^ > 1.6) 
were fitted with rmodei — 1 + C by varying parameters 
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FIG. 4: Left: pair-density net ratio rmodei[X{pti), X{pt2)] — 
1 for model fit to mid-central (b) Au-Au collisions. Right; 
residuals (data — model) for mid-central collisions. 



ris, "A and C (offset). Parameters /3o = 5 GeV^^ and 
Too — ttLt^ were fixed by the (pion dominated) inclusive 
single-particle pt spectrum forpt < 1 GeV/c. The fits are 
insensitive to the absolute value of l/rifiuct'i its value was 
fixed as follows. Parameter 1/n when fitted to the sin- 
gle particle rrit spectrum [l^ , using an analog of Eq. (j4|) 
with Ufiuct replaced by n, accounts for the deviation be- 
tween the measured distribution and e"''™*. In general, 
both collective radial expansion velocity [11] and effec- 
tive temperature fluctuations contribute to the curva- 
ture (decreasing slope) of the rrit spectrum relative to 
Boltzmann reference e"''™' at increasing rrit shown by 
the dashed curve in Fig.[3l^c). Both contributions are in- 
cluded in parameter n in Eq.(4), when fitted to the single 
particle distribution, resulting in an apparent variance, 
1 /n, given by an incoherent sum of contributions from ra- 
dial fiow, l/ufiow, and effective temperature fluctations, 
'^/rifiucu where 1/n = l/ufiow + '^/rifiuct- However, for 
the effective temperature fluctuation model developed in 
the preceding section only component l/ufiuct is rele- 
vant to the 2D Levy saddle flt but it is not accessible 
because fits to correlation data (r — 1) poorly constrain 
absolute quantities l/n^ and I/tt-a. However, differences 
A(l/n)5],A are well determined by the saddle curvatures, 
nearly independently of the assumed value of l/rifiuct 
in r model- The maximum value for l/rifiuct corresponds 
to 1/n = 1/13 in the no-fiow limit, l/rifiow = 0, where 
n = 13 is obtained from the Levy distribution fit to the 
single particle rrit spectrum [T8| . The minimum value of 
0.0009 corresponds to that necessitated by the fitted val- 
ues of A(l/n)A in Table [J in the limit 1/nA ^ 0. The 
fits were insensitive to variations of l/nj;„ct in this range, 
intermediate value l/njiuct — 0.03 near the center of the 
allowed range provided stable A(l/n)s and A(l/n)A fit 
values. Best- fit parameters and x^/DoF for the saddle 
fits are listed in Table HI The model function and resid- 
uals for the fit to centrality (b) are shown in Fig. [D 

Two-dimensional saddle-fit residuals, as in Fig. H (right 
panel), are approximately constant along directions par- 
allel to the X{pt)A — X{Pti) — X{pt2) axis for each value 
of X{pt)^ and are small for X{pt)j: < 1.5. The Levy 



TABLE 1: Parameters and fitting errors (only) for 2D veloc- 
ity/temperature fluctuation model for each centrality bin, (a) 
- (d) (central - peripheral as in Fig. [ij . Errors (last column) 
represent fitting uncertainties. Systematic errors are 7-12% 
[301 . Mean multiplicities of used particles in the acceptance, 
N, are listed for each centrality bin. Quantities S (last row) 
are correction factors for contamination and tracking ineffi- 
ciency 



centrality 


(d) 


(c) 


(b) 


(a) 


error" (%) 


N 


115.5 


424.9 


790.2 


983.0 




C X lO"* 
A(l/n)E X 10"* 
A(l/n)A X 10'* 
A{l/n)tot X 10^ 

xVdof 


-11.6 
3.54 
-8.61 
12.2 

348 
286 


-0.820 
0.611 
-3.33 
3.95 

313 

286 


0.787 
0.183 
-2.53 
2.71 

475 
286 


0.750 
0.118 
-2.04 
2.16 

402 
286 


6-14 
6-24 
6-3 


s 


1.19 


1.22 


1.25 


1.27 


8^ 



"Range of fitting errors in percent from peripheral to central. 
''Systematic error. 



temperature fluctuation model adequately describes the 
saddle structure. Residuals from the fit for mid-central 
events (b) are shown in Fig. [5] (left panel) projected onto 
sum variable X{pt)s- Errors are included in the data 
symbols and are smaller than those in Fig. [4] (right panel) 
due to bin averaging. Residuals for other centralities are 
similar, but differ in amplitude. We hypothesize that this 
residual structure is due to correlated final-state hadrons 
associated with initial-state semi-hard parton scatter- 
ing M- 

Centrality dependences of efficiency-corrected model 
parameters 5iVA(l/n) [4l[, which determine saddle- 
shape correlation amplitudes in Fig. [1] are shown in 
Fig. O (right panel). The linear trends suggested by 
the solid lines are notable. Multiplication by factor 
SN estimates correlation amplitudes per final state pri- 
mary particle as discussed below. Centrality measure 
v estimates the mean participant path length as the 
average number of encountered nucleons per partici- 
pant nucleon in the incident nucleus. For this anal- 
ysis 1^ EE 5.5{N/Noy/^ ~ 5.5 {Npart/Npart,^axy^^ ^ 

'^Nbin/Npart, based on Glauber- model simulations where 
Npart (Nbin) is thc numbcr of participant nucleons (bi- 
nary collisions). 

The reasons for multiplying the parameters in Ta- 
ble H] by iS and N are the following. Multiplication of 
(f — 1) by iV yields the density of correlated pairs per 
final-state particle plj . typically 0(1) for all centralities. 
iV(f — 1) would be independent of centrality if Au-Au 
collisions were linear superpositions of p-p collisions (par- 
ticipant scaling) because the amplitude of the numerator 
of (f — 1), which is proportional to thc density of corre- 
lated pairs, would scale with participant number, or in 
this model with N, while the denominator is proportional 
to 7V^. Therefore variation of N{f — 1) with centrality 
directly displays the effects of those aspects of Au-Au 
collisions which do not follow naive p-p superposition. 
Factor S is defined as the ratio of true, primary particle 
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FIG. 5: Left: Residuals from 2D Levy saddle fit to mid-central 
(b) data in Fig. [1] projected onto sum variable X{pt)s — 
X{pti) + X{pt2)- Right: Efficiency-corrected per-particle 
saddle-curvature measures on centrality v. 5iVA(l/n)E 
(dots), —SN A{l/n)A (triangles) and SN A{l/n)tot (open 
circles). Data symbols include fitting errors only Solid 
lines are linear fits. 



yield {i.e., 100% tracking efficiency and no background 
contamination) estimated for these data in Ref. [ISj di- 
vided by the actual multiplicity used in this analysis cor- 
rected for the ~ 7% background contamination. S is 
essentially the reciprocal of the charged-particle tracking 
efficiency, specific for the present analysis. Multiplica- 
tion by factor SN of the parameters in Table [T] there- 
fore estimates the correlation amplitudes per final-state 
particle for 100% tracking efficiency and no background 
contamination, assuming the measured correlations in- 
clude background-primary particle correlations half-way 
between the limits described in Sec. IIII CI The uncer- 
tainty in extrapolating to the true primary particle yield 
is estimated to be 8%, most of which is due to the 7% 
systematic uncertainty in the measured charged hadron 
yield The combined systematic uncertainty for the 
efficiency corrected amplitudes is from ±11 — 14% across 
the X{pti) vs X{pt2) space. 



VI. DISCUSSION 

Correlations onpt have two main components, a saddle 
shape and a peak at higher pt. By measuring the saddle 
curvatures we infer the relative covariance of two-point 
distribution g2{(3i,l32) and hence the average two-point 
correlation amplitude of the temperature/velocity struc- 
ture of the composite particle source. We now consider 
possible dynamical origins of that structure. 

The analysis of the saddle-shape produces accurate re- 
sults for relative variance differences A(l/n)s = (cr^^, — 
(T^)//3g, A(l/n)A = (o-^^ - cr|)//3o, and the correspond- 
ing A(l/n)tot = (cr^j, — '^^j^)/0o fo'^ effective tempera- 
ture fluctuations. The measurements do not constrain 
the absolute magnitudes of the individual variances, cr|^ 
and a^^ . The minimum possible values, consistent 
with the saddle-shape conditions and the single-particle 
rrit spectra, correspond to ct^^ = and l/rifiuct — 



-A(l/n,)A, resulting in (T/3//3o = ctt/Tq = 1.4% to 
2.9% global event-to-event temperature/velocity fluctua- 
tion from central to peripheral collisions, respectively. In 
this case l/n = l/ufiow and global temperature/ velocity 
fluctuations contribute negligibly to the upward curva- 
ture of the dN/nitdmt spectrum. The maximum values 
for the variances correspond to l/ufiow = 0, resulting 
in l/njiuct = ^/n and ap/Po = or/To = i/l/n = 30%, 
where CT/j^ ~ cr^g, corresponding to 30% local temper- 
ature/velocity fluctuations within each event, a signifi- 
cantly non-equilibrated system. Thus, local temperature 
variation could range between and 30%. One can ask 
what is the source of the fluctuating effective tempera- 
ture, and is local source velocity rather than temperature 
a more appropriate quantity? 

Given the correlation peaks at higher pt it is reason- 
able to offer the hypothesis that the saddle-shape corre- 
lation structure in Fig. [T] results from in-medium modi- 
fication, specifically momentum dissipation on (j>ti,Pt2) 
of a two-particle distribution from fragmenting, semi- 
hard scattered partons in the initial-stage of the colli- 
sion. Since no selection was made on leading particle or 
high-pt "trigger" particle for these data we refer to the 
hadrons associated with a semi-hard, initial-state scat- 
tered parton as a minijet (ol. |4^. Minijet production in 
Au-Au collisions should increase approximately linearly 
with Nbin ^M, H] while the subsequent momentum dissi- 
pation should monotonically increase with greater mini- 
jet production. Correlation amplitudes per final state 
particle (the latter approximately proportional to Npart) 
should therefore increase monotonically with mean par- 
ticipant path length u = Nhm/ (Npart/'i), thus providing 
a basis for experimental tests of this hypothesis. 

The linear trends in Fig. [5] (right panel) therefore sup- 
port, but do not require, a minijet — momentum dissipa- 
tion mechanism for the observed correlations on pt. In 
Fig. [5] we also observe 1) reduced curvature along the 
sum direction and 2) increased curvature along the dif- 
ference direction which may represent respectively trans- 
port of semi-hard parton structure to lower pt and a more 
correlated bulk medium. Given a minijet interpretation 
of S N A{1/ n)tot, the combined trends 1) and 2) repre- 
sent strong evidence for increased parton dissipation in 
the more central Au-Au collisions. The present results 
complement the observed suppression of high-pt spectra 
{Raa) [H, Hi and suppression of large angle trigger- 
particle - associated-particle conditional distributions on 
A(j) [1, § in central Au-Au collisions at RHIC. It is very 
likely that the lower-pt fluctuations and correlations re- 
ported here are, at least in large part, a consequence of 
the processes which lead to the above suppressions at 
higher-pt . 

It is important to note that these correlations on trans- 
verse momentum observed at relatively low pt reveal 
nominally 'soft' structure in relativistic heavy ion col- 
lisions which scales with the number of binary collisions 
Nhin, whereas a low-pt inclusive quantity such as mul- 
tiplicity scales with participant number Npart- Binary- 
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collision scaling is conventionally thought to be an aspect 
of high-pt physics and initial-state scattering. This analy- 
sis suggests that substantial effects of initial-state parton 
scattering are manifest at low pt in more central heavy 
ion collisions. 



VII. SUMMARY 

In conclusion, the dynamical origins of excess (pt) 
fluctuations in Au-Au collisions at RHIC are studied 
in this analysis of two-particle correlations on {pti,pt2)- 
The velocity/temperature structure of heavy ion colli- 
sions suggested by these correlations is unanticipated by 
theoretical models [l^, [ll|. Lacking in these mod- 
els is the simultaneous inclusion of hard scattering in 
the initial state with subsequent medium modification of 
the fragmentation function and/or interactions between 
the medium and the hadrons associated with the scat- 
tered partons. Nevertheless it seems plausible to inter- 
pret the observed correlations on {pti,Pt2) as resulting 
from this sort of semi-hard parton scattering and subse- 
quent medium modified fragmentation and/or associated 



hadron distributions on pt in the more central Au-Au 
collisions. In this picture, with increasing centrality the 
transverse momentum associated with the two-particle 
fragment distribution from initial-state semi-hard parton 
scattering is shifted to lower pt , asymptotically approach- 
ing a form consistent with random velocity/temperature 
variations (Levy saddle) as a manifestation of substan- 
tial but incomplete equilibration. These newly-observed 
Pt correlations may thus reveal minijet dissipation in the 
medium produced by Au-Au collisions at RHIC. 
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